A geometric generalization of field theory to manifolds of arbitrary dimension 



Kay Jorg Wiese 
Fachbereich Physik, Universitat GH Essen, 45117 Essen, Germany 

Mehran Kardar 

Department of Physics, MIT, Cambridge, Massachusetts 02139, USA 

We introduce a generalization of the O(N) field theory to iV-colored membranes of arbitrary 
inner dimension D. The O(N) model is obtained for D 1, while N — *• leads to self-avoiding 
tethered membranes (as the O(N) model reduces to self-avoiding polymers). The model is studied 
perturbatively by a 1-loop renormalization group analysis, and exactly as N — + oo. Freedom to 
choose the expansion point D, leads to precise estimates of critical exponents of the O(N) model. 
Insights gained from this generalization include a conjecture on the nature of droplets dominating 
the 3d-Ising model at criticality; and the fixed point governing the random bond Ising model. 
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Field theories have strong connections to geometrical 
problems involving fluctuating lines. For example, sum- 
ming over all world-lines representing the motion of par- 
ticles in space-time, is the Feynman path integral ap- 
proach to calculating transition probabilities, which can 
also be obtained from a quantum field theory. Another 
example is the high-temperature expansion of the Ising 
model, where the energy-energy correlation function is a 
sum over all self-avoiding closed loops which pass through 
two given points. Generalizing from the Ising model to N 
component spins, the partition function of a correspond- 
ing O(N) "loop model" is obtained by summing over all 
configurations of a gas of closed loops, where each loop 
comes in N colors, or has a fugacity of N. In the limit 
N — > 0, only a single loop contributes, giving the parti- 
tion function of a closed self-avoiding polymer 

There are several approaches to generalizing fluctuat- 
ing lines to entities of other internal dimensions D. The 
most prominent examples are string theories and lattice 
gauge theories, which both describe D = 2 world sheets 
[]2j. The low temperature expansion of the Ising model 
in d dimensions also results in a sum over surfaces that 
are d—1 dimensional. Each of these generalizations has 
its own strengths, and offers new insights on field theory. 
Here we introduce a generalization based on a class of 
-D-dimensional manifolds called "tethered" (or polymer- 
ized) membranes, which have fixed internal connectivity, 
and are the simplest generalization of linear polymers || . 

Tethered manifolds are a quite successful generaliza- 
tion of self-avoiding polymers embedded in d-dimensional 
external space: Simple power counting indicates that 
self- avoidance is relevant only for dimensions d < d c — 
AD I '(2 - D), making possible an e = 2D - d(2 - D)/2 ~ 
(d c (D) — G?)-expansion, which was first carried to 1-loop 
order in Refs. Q]. Following more rigorous analysis of this 
novel perturbation series g , recently 2-loop calculations 
were performed for 1 < D < 2 ||. To obtain results for 
polymers or membranes, one has the freedom to expand 
about any inner dimension D, and the corresponding up- 
per critical dimension of the embedding space 0. This 
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FIG. 1. Schematic depiction of the model, and its limits. 

freedom can be used to optimize calculation of critical ex- 
ponents, and when applied to polymers, the results are 
clearly better than those from the standard expansion. 

Our basic idea is to generalize the high temperature ex- 
pansion of the 0{N) model, from a gas of self-avoiding 
loops of fugacity N, to a similar gas of closed fluctu- 
ating manifolds of internal dimension D. The primary 
goal is to obtain a novel analytical handle on the field 
theory for D — 1, and we do not insist that the mod- 
els for general D correspond to any physical problem. 
Given this caveat, the generalization is not unique. En- 
couraged by its success in polymer theory, we study teth- 
ered manifolds, and in addition restrict ourselves to the 
genus of hyper-spheres. We have chosen hyperspheres as 
they have no additional anomalous correction exactly at 
D = 2. The resulting manifold theory depends on two 
parameters N and D, whose limiting behaviors reduce to 
well known models, as indicated in Fig. |l| 

We construct the M(N,D) manifold model perturba- 
tively, starting with the well known high-temperature ex- 
pansion (loop-model) of the O(N) field theory as 
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The first term in this sum is the partition function of a 
non-interacting polymer loop, now generalized to a D- 
dimensional manifold, fixed at the origin in d dimensions 
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The integral runs over all sizes f2, weighted by a chemical 
potential t, while the factor 



S d (r(x )) 
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is the probability that the membrane is attached to a 
given point in space ||. The above generalization de- 
pends on a function c(D) , related to the relative strengths 
of self-avoidance between parts of the same manifold, and 
between different manifolds. Any choice of c(D) which 
reproduces the polymer partition function with c(l) = 1, 
is acceptable. In the remainder, we will mainly focus on 
c(D) = D, equivalent to the integral dO over all scales. 

Subsequent terms in Eq. (Q) correct for the intersec- 
tions (symbolized by a dashed line) between manifolds. 
Configurations involving intersections between different 
manifolds, or self- intersections of the same manifold (e.g., 
the second and third terms in the above expansion, re- 
spectively) have to be subtracted from the partition func- 
tion. The contributions from these terms result in di- 
vergencies which have to be removed by renormalization 
(for details see Ref. M). The final result is the 1-loop 
expansion of the exponent v(D, N, d) (for the divergence 
of the correlation length at a critical point), given by 
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This expression reduces to the well-known e-expansion 
H around d = 4 for lines (D — 1), while the N — > limit 
reproduces the result for self-avoiding manifolds 0] . 

The perturbative series can also be summed exactly in 
the large N limit, where the ambiguity associated with 
c(D) is removed B, to give the exponent 
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This generalizes the well-known result of v = l/(d — 2) 
for the 0(N) model §. 

To demonstrate the utility of this generalization, we 
next estimate the exponent v in physical dimensions, by 
expanding about a point (Do, do) on the critical curve 
e(D , d ) — 0. The simplest scheme is to extrapolate to- 
wards the physical theories for D — 1, 2 and d = 2, 3, • • •, 
with the expansion parameters D — Dq and d — do ■ How- 
ever, this choice is not optimal (see Ref. |(|), and better 
results are obtained with expansion parameters D c (d) = 
2d/(4 + d) and e(D, d) = 2D - d{2 - D)/2d. Further- 
more, it is better to expand vd or v(d + 2) rather than v. 
The first is an expansion of the form vd = 2D + a(D)e, 
about the mean-field result vmf — 2D/d, also known 
as Gaussian variational approximation pQ| . The sec- 
ond is a similar expansion about the Flory expression 
^Flory = (2 + D)/(2 + d). 
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FIG. 2. Extrapolations for the exponent v of the O(N) 
model in d — 3, using the expansion of vd with c(D) — D. The 
dashed lines represent the best known values from Ref. ||. 

After selecting one of these schemes, Eq. (|]) is re- 
expressed in terms of the chosen expansion parameters. 
However, we are still free to choose the expansion point 
along the critical curve, which then fixes Dq. As Dq 
is varied, different values for v are obtained, as ploted 
in Fig. |^. The criterion for selecting a particular value 
from such curves is that of minimal sensitivity to the ex- 
pansion point, and we thus evaluate v at the extrema. 
The broadness of the extremum then provides a mea- 
sure for the goodness of the result, and the expansion 
scheme. Although we examined several such curves, 
only a selection is reproduced in Fig. ||. Our results 
are clearly better than the standard 1-loop expansion of 
v = l/2 + (JV + 2)/[4(JV + 8)]. 

While providing good estimates for exponents of the 
O(N) model is certainly a benefit, the generalization to 
M(N,D) should offer insights beyond the standard field 
theory. Also, for the approach to be more generally ap- 
plicable, we should show that similar generalizations are 
possible for other field theories. In the rest of this article 
we shall demonstrate that these goals are indeed feasible. 

Our starting point was the high temperature expan- 
sion of the O(N) spin model, which naturally leads to 
a sum over iV-colored loops (D = 1); motivating the 
later generalization to M(N,D). For the Ising model 
(N = 1), a different geometrical description is obtained 
from a low temperature expansion: Excitations to the 
uniform (up or down pointing) ground state are droplets 
of spins of opposite sign. The energy cost of each droplet 
is proportional to its boundary, i.e. again weighted by a 
Boltzmann factor e~ . Thus, a low temperature series 
for the d-dimensional Ising partition function is obtained 
by summing over closed surfaces of dimension D = d — 1 . 
For d — 2, the high and low temperature series are sim- 
ilar, due to self-duality. For d = 3, the low temperature 
description is a sum over surfaces. 
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FIG. 3. Test of Eq. (|) for Ising models in d = 2.5 and 
3. The upper curves are from the high temperature repre- 
sentation (D = 1), while the lower curves are from the low 
temperature description (D = d — 1). The exponent v is esti- 
mated from the maximum of each curve, which are obtained 
by extrapolating vd with c(D) — D, linearizing in N, and 
dividing the result by dD. 

The non-trivial question is regarding the types of sur- 
faces which dominate the above sum. Since there is no 
constraint on the internal metric, it may be appropriate 
to examine fluid membranes. However, there is currently 
no practical scheme for treating interacting fluid mem- 
branes, and the excluded volume interactions are cer- 
tainly essential to the problem. It is known that con- 
figurations of a single surface for N — 0, self-avoiding 
or not, are dominated by tubular shapes (spikes) which 
have very large entropy pH] . Such "branched polymer" 
configurations are very different from tethered surfaces. 
However, for N ^ 0, it may be entropically advantageous 
to break up a singular spike into a string of many bubbles. 
If so, describing the collection of bubbles by fluctuating 
hyper-spherical (tethered) manifolds may not be too off 
the mark [|l6|. We can test the validity of this conjec- 
ture by comparing the predictions of the dual high and 
low-temperature descriptions. 

Singularities of the partition function are characterized 
by the critical exponent a{D,d, TV), or (using hyperscal- 
ing) through In Z s i ng ui ar ~ \t— t c \ vd ' D . The equality of the 
singularities on approaching the critical point from low 
or high temperature sides, leads to a putative identity 

v{l,d,l) = -^—v{d-l,d,l) ■ (6) 
d — 1 

Numerical tests of the conjecture in Eq. (||) are presented 
in Fig. ||. The extrapolated exponents (the maxima of the 
curves) from the dual expansions are in excellent agree- 
ment. Nevertheless, higher-loop calculations would be 
useful to check this surprising hypothesis. 

The simplest extension of the 0(N) model breaks the 
rotational symmetry by inclusion of cubic anisotropy [fL2| . 
In the field theory language, cubic anisotropy is repre- 
sented by a term uJ2i s ti m addition to the usual in- 
teraction of by] - sfs'j. In the geometric prescription of 
high temperature expansions, the anisotropic coupling u 
acts only between membranes of the same color, while 
the interaction b acts irrespective of color. 
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FIG. 4. Regions with different RG-fiow patterns in the 
(N, D)-plane (top), and the corresponding RG- flows (bot- 
tom); shaded regions are unstable. 



Stability of the system of colored membranes places 
constraints on possible values of b and u. To avoid col- 
lapse of the system, energetic considerations imply that 
if u < 0, the condition u + b > must hold, while if 
u > 0, we must have u + Nb > JT^J§]. The above sta- 
bility arguments are expected to be modified upon the 
inclusion of fluctuations. A well-known example is the 
Coleman- Weinberg mechanism fll^| , where the RG flows 
take an apparently stable combination of b and u into an 
unstable regime, indicating that fluctuations destabilize 
the system. In the flow diagrams described below, we also 
find the reverse behavior in which an apparently unsta- 
ble combination of b and u flows to a stable fixed point. 
We interpret this as indicating that fluctuations stabilize 
the model, a reverse Coleman- Weinberg effect, which to 
our knowledge is new. We have shaded in grey, the un- 
physical regions in the flow diagrams of Fig. ^. As in 
their O(N) counterpart, the RG equations admit 4 fixed 
points: The Gaussian fixed point with b* G = u* G = 0; 
the Heisenberg fixed point located at b* H ^ 0, u* H = 0; 
the Ising fixed point with b*j = 0, u*j =/= 0; and the cubic 
fixed point at b* c ^ 0, u* ^ 0. Furthermore, as depicted 
in Fig. [|, there are six different possible flow patterns. 
In the O(N) model, the flows in (i) and (ii) occur for 
N < 4 and N > 4, respectively. The other four patterns 
do not appear in the standard field theory, as is apparent 
from their domain of applicability in the (TV, D) plane on 
Fig. ^. Note that there are two stable fixed points in three 
out of these four cases. 




FIG. 5. Extrapolations of v from the expansion of vd with 
c(D) = D, for the O(N) model in d = 3. Exponents at the 
Heisenberg fixed point for TV = are compared to those of the 
Ising and cubic fixed points. The crossing of the latter curves 
yields an estimate of v = 0.6315 for the 3d Ising model. 

The N — > limit of the above models is interesting, 
not only because of its relevance to self-avoiding poly- 
mers and membranes, but also for its relation to the Ising 
model with bond disorder. The latter connection can be 
shown by starting with the field theory description of the 
random bond Ising model, replicating it N times, and av- 
eraging over disorder [u| . The replicated system is con- 
trolled by a Hamiltonian with positive cubic anisotropy 
u, but negative b = —a (a is related to the variance of 
bond disorder). From the "Harris criterion" |l4[ , new 
critical behavior is expected for the random bond Ising 
system. But in the usual field theory treatments [ p"3| , 
there is no fixed point at the 1-loop order. In our gener- 
alized model, this is just the borderline between cases (i) 
and (iii). However, we now have the option of searching 
for a stable fixed point by expanding about any D =/= 1. 
Indeed, for N = and 1 < D < 1.29, the cubic fixed 
point lies in the upper left sector (u > and b < 0) and 
is completely stable, as in flow pattern (iii). 

The extrapolation for v at the cubic fixed point is plot- 
ted in Fig. ||, where it is compared to the results for the 
Heisenberg and Ising fixed points. The divergence of v 
on approaching D — 1 from above, is due to the cubic 
fixed point going to infinity as mentioned earlier. Upon 
increasing D, the Ising and cubic fixed points approach, 
and merge for D = 1.29. For larger values of D, the cu- 
bic fixed point is to the right of the Ising one (b* > 0), 
and only the latter is stable. Given this structure, there 
is no plateau for a numerical estimate of the random 
bond exponent z^do, and we can only posit the inequal- 



ity ^do > Rising- While this is derived at 1-loop order, it 
should also hold at higher orders since it merely depends 
on the general structure of the RG flows. One may com- 
pare this to four loop calculations of the random bond 
Ising model Jig] , which are consistent with a = 0, i.e. at 
the border-line of the Harris criterion [0, with v = 2/3. 
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